Resonance tongues and patterns in periodically forced reaction-diffusion systems 
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Various resonant and near-resonant patterns form in a ligiit-sensitive Belousov-Ziiabotinsky (BZ) reaction 
in response to a spatially-homogeneous time-periodic perturbation with light. The regions (tongues) in the 
forcing frequency and forcing amplitude parameter plane where resonant patterns form are identified through 
analysis of the temporal response of the patterns. Resonant and near-resonant responses are distinguished. 
The unforced BZ reaction shows both spatially-uniform oscillations and rotating spiral waves, while the forced 
system shows patterns such as standing-wave labyrinths and rotating spiral waves. The patterns depend on 
the amplitude and frequency of the perturbation, and also on whether the system responds to the forcing near 
the uniform oscillation frequency or the spiral wave frequency. Numerical simulations of a forced FitzHugh- 
Nagumo reaction-diffusion model show both resonant and near-resonant patterns similar to the BZ chemical 
system. 

PACS numbers: 82.40.-g, 05.45Xt, 05.65+b, 05.45.-a 



I. INTRODUCTION 

An oscillator forced by a periodic external perturbation en- 
trains to the forcing for certain values of the perturbation fre- 
quency and amplitude. This behavior is observed in a wide 
range of biological, chemical and physical systems, for exam- 
ple, in circadian rhythms such as the sleep-wake cycle forced 
by the sun 1 1 ], in the tips of chemical spiral waves forced with 
light 4^^^, and in arrays of Josephson's junctions [TJ. 

The entrainment to the forcing can take place even when 
the oscillator is detuned from an exact resonance 
In this case, a periodic force with a frequency // shifts the 
oscillator from its natural frequency, /o, to a new frequency, 
fr, such that // /fr is a rational number m:n. When the forcing 
amplitude is too weak this frequency adjustment or locking 
does not occur; the ratio // / fr is irrational and the oscillations 
are quasi-periodic. In dissipative systems frequency locking is 
the major signature of resonant response. Nearly conservative 
systems show in addition a large increase in the amplitude of 
oscillations. 

The response of a two-dimensional array of coupled non- 
linear oscillators, or of a two-dimensional oscillating field is 
much less well understood. For a periodically forced single 
oscillator, the structure in the parameter plane of the forcing 
frequency // and amplitude / contains many universal fea- 
tures identified with frequency locking 1 11], but it is unclear 
if these features persist or change in the more complex case 
of spatially extended systems where patterns can form. Pat- 
terns resulting from time-periodic forcing [1_2, 13^ 1 4, _15. 16. 
[13 El [HIS Elil, spatially periodic forcing ElTand 
global feedback fz?, '25^,^0] have been studied in the past, but 
a whole phase diagram showing a resonance structure for spa- 
tially extended systems has not been reported. Moreover, it is 
not even clear to what extent the familiar concept of resonance 



applies to spatially extended systems. 

To investigate a phase diagram showing multiple reso- 
nances in spatially extended systems, we applied periodic 
perturbations to a light sensitive form of the Belousov- 
Zhabotinsky (BZ) chemical reaction-diffusion system and 
measured the temporal response and pattern formation. In the 
course of this investigation a new mode of response has been 
identified. Patterns may fail to lock to the forcing frequency 
but still respond by showing an m-peaked distribution of the 
oscillation-phase as in resonant patterns. We refer to this re- 
sponse mode as "near-resonant." 

In addition to the nonuniform distribution of the oscilla- 
tion phase, resonant and near-resonant, patterns can also be 
characterized by the shape of the phase in the complex phase 
plane. The phase of unforced spirals has a circular shape in 
the complex phase plane but forcing breaks the circular sym- 
metry. At high enough forcing this is visible as a m-fold sym- 
metry in the phase plane. Examples of m:l patterns observed 
in the BZ system, where m:l=////o and m 2,3,4 are shown 
in Fig[2 Unlike the single oscillator case, a spatially extended 
system can exhibit phase waves and other phase patterns. In 
Fig.^each pattern is shown in two representations: in the real 
space x — y plane, and in the complex phase plane. 

In this paper we construct an experimental phase diagram 
in the forcing frequency and amplitude parameter plane of 
resonant and nearly resonant m:n responses and identify the 
pattern types that lead to the two responses. The experi- 
mental set-up and determination of resonance tongues are 
described in Sec. [Il] and Sec. [IllJ respectively. The qualita- 
tively different patterns observed in the experiments are pre- 
sented in Sec. IIVI A forced reaction-diffusion model, mod- 
ified FitzHugh-Nagumo equations, is introduced in Sec. [V] 
which is followed by a discussion of the model and the in- 
sight it provides into the mechanisms of pattern formation in 




FIG. 1: (Color online) Patterns in the periodically forced BZ reaction, (a)-(d) show a 11.5 x 11.5 mm size region of patterns formed at 
different forcing frequencies and amplitudes. The data are processed to show the pattern only near the sub-harmonic response frequency. 
Patterns (a) and (c) are frequency-locked (resonant). Patterns (b) and (d) are near-resonant but not frequency-locked. The frames (e)-(h) show 
the same data represented in the complex phase plane, see 1 2 1 ] . The colors in the phase plane vary with the angle [0, 2n] and are the key to the 
pattern images above. See the text for details. (a,e) 2:1 two phase standing-wave pattern; (b,f) 3:1 three phase rotating spiral (c,g) 3:1 three 
phase standing-wave pattern; (d,h) 4:1 four phase rotating spiral. Chemical conditions are given in Ii27i1 . 



the experimental system. A general discussion and summary 
of the results are given in Sec. lVll 



II. BELOUSOV ZHABOTINSKY CHEMICAL SYSTEM 

The oscillatory chemical reaction occurs in a thin, porous- 
glass membrane (0.4 mm thick, 25 mm in diameter), which 
is in contact on each side with continuously refreshed reser- 
voirs of reagents for the ruthenium-catalyzed Belousov- 
Zhabotinksy (BZ) reaction |28]. Each reservoir is well-stirred 
and the reagents diffuse from them into the membrane where 
they react. The chemical concentrations in the reservoirs are 
given in Ref. f27','2^. Visualization of the patterns is achieved 
using a low intensity tungsten lamp, which measures the op- 
tical density of the concentration patterns in the membrane 
without affecting the chemical reactions. 

For the chemical concentrations used in the present ex- 
periments the unforced system exhibits rotating spiral pat- 
terns. The patterns are sustained indefinitely in time because 
the reaction products leave the membrane by diffusion into 
the reservoirs, and reservoir concentrations are maintained by 
continuous feeds. We used two different sets of chemical con- 
ditions ll27ll29il . one creating spirals with a higher frequency 
ifs = 0.072 Hz) and one creating spirals of a lower frequency 
(fs — 0.020 Hz). Examples of the unforced spiral waves for 
both sets of chemical conditions are shown in Fig.|2] 

In addition to the spiral frequency, the BZ system has an- 
other natural frequency: the unforced spatially-homogeneous 
oscillation frequency. Since perturbations always lead ulti- 
mately to the formation of spiral waves in the membrane, we 




FIG. 2: Unforced spiral patterns for the two sets of chemical con- 
ditions used for the experiments presented, (a) Spiral waves with a 
shorter period, [BrMA]=0.220 M, [firO3"]=0.230 M |27]. (b) Spi- 
ral waves with a longer period, [BrMA]=0300 M, [SrO3"]=0.136 
M 1 29]. The images show a 9 mm x 9 mm region of the pattern. Dark 
regions correspond to high concentrations of Ru(II). Taken from 1 12]. 
All data reported in this paper are taken under the conditions for the 
spirals on the left, except for Fig.|4lb). 



determined the homogeneous oscillation frequency /o by the 
following method. The membrane was exposed to a spatially 
uniform high-intensity pulse of light for 30 s, which resets the 
system so the entire membrane is oscillating with the same 
frequency and phase. The frequency of this spatially uniform 
oscillation, determined by the chemical kinetics, was found in 
our previous work to be essentially independent of the chemi- 
cal concentrations used in our study, /o — 0.020 Hz 1 12]. The 
uniform oscillations eventually evolve to rotating spiral waves 
which fill the system. However, the frequency of the spiral 
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waves does depend on the chemical conditions. The ramifica- 
tions of this dependence will be discussed in the next Section. 

m. MULTIPLE RESONANCE TONGUES 

The chemical reaction is forced by illuminating it with spa- 
tially homogeneous light that is periodically blocked 
I21I1 ; the durations of the illuminated and blocked portions of 
each cycle are equal, i.e., the intensity modulation is a square 
wave. To investigate the temporal response of the patterns 
to the forcing we varied the intensity / and frequency of the 
periodic light forcing //. The light intensity / is the control 
parameter for the forcing amplitude. We examined the tem- 
poral response to determine which, if any, tongue the pattern 
belonged to, and examined the existence, shape, and ordering 
of the resonance tongues (see Sec.lllH. 

A. Determining Temporal Resonance 

The experimental data were collected as a time-sequence of 
pattern snapshots. The natural oscillation period of the reac- 
tion for the conditions used was about 50 s. Typically, images 
were recorded every 2 s for one hour and a central 240 x 240 
pixel (23 mm x 23 mm) region of the pattern was analyzed. 
The Fourier transform of the time series for each pixel was 
calculated to obtain an average power spectrum for the pat- 
tern. Figure|3]shows a typical averaged power spectrum for a 
resonant pattern when the forcing frequency // was twice the 
uniform oscillation frequency /q. The largest subharmonic 
frequency peak appears at / = ff/2, as indicated by the verti- 
cal line. 



B. Tongues 

Using the method described in the previous section, we 
found tongues in the forcing parameter space where the pat- 
tern responds at or near m:n resonances. We obtaine d a p hase 
diagram for each of the two chemical conditions [23 l2^ . 
shown in Fig. EIa),(b). If the peak of the strongest mode 
subharmonic to the forcing was within ±3% of the forcing 
frequency, we considered the pattern to be responding to the 
forcing, and it is included within an m:n resonance region in 
Fig. 13 This criterion is consistent with the observation of 
m-peaked distributions of the oscillation phase. Some of the 
patterns meeting this criterion are resonant while others are 
near-resonant, i.e., they are quasi-periodic patterns with an m- 
peaked phase distribution. 

We varied the forcing frequency and intensity in the ex- 
periments and explored the temporal resonant response as we 
moved through the parameter space, and the results are shown 
in Fig. 0] Each symbol type represents a different m:n res- 
onance. The curves in Fig. 0] are drawn to guide the eye to 
the tongues in the ff-I plane with different m:n responses. 
Only the largest resonance tongues are plotted. In addition 
to the m:\ tongues (and the 4:3 and 3:2 tongues) shown, we 
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FIG. 3: Averaged power spectrum for a 2:1 resonant labyrinthine 
pattern. The response frequency has been normalized to the forcing 
frequency ff. The peak at 0.5 is the subharmonic response to the 
forcing. The large peak at the forcing frequency (/ = ff) is due to 
our imaging method, which captures part of the forcing light output. 
In this paper we only consider responses subharmonic to the forcing 
frequency (//// < 1)). 



observed several higher order m:n states (e.g., 5:7, 5:1, 6:1, 
10:1), which spanned control parameter ranges too narrow to 
be maintained. In all cases the different m:n tongues were or- 
dered in a Farey sequence, similar to the Devil's staircase or- 
dering of resonance tongues for two coupled oscillators fsoll 
and for the homogeneous BZ reaction |31, 32]. 

We investigated two different chemical conditions. The 
chemical conditions that yield 0.072 Hz spirals have m:n 
tongues that bend toward higher frequency as the light inten- 
sity / is decreased [Fig. Ua)]. For the chemical conditions 
that yield 0.020 Hz spirals, the tongues do not bend much 
at low frequency [Fig.^b)]. The bending of the tongues is 
caused by a shifting from /o = 0.02 Hz resonance at high forc- 
ing intensity (the uniform oscillation frequency) to the near- 
resonant response of the spiral wave frequency fs for lower 
forcing intensity. Since the spiral wave frequency for the data 
in Fig. |4lb) is the same as the uniform oscillation frequency, 
the tongues do not bend in that case lll2ll33ll . 



C. A quantitative measure of patterns 

Spatial Fourier transforms and correlation functions do not 
capture the temporal aspects of the patterns and did not dif- 
ferentiate the data well because the patterns were often com- 
prised of multiple wavelengths and orientations. Therefore, 
instead of computing spatial Fourier transforms, we analyzed 
the temporal Fourier transform calculated for each point in 
the pattern O 111 El]. The power spectrum of the signal, 
averaged over the spatial pattern, gives information about the 
strongest frequency response. Additionally, the data were fil- 
tered to keep the strongest response and then inverse Fourier 
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FIG. 4: (Color online) The largest m:n tongues observed in the 
frequency-intensity plane of the spatially extended BZ system. 
Chemical conditions are those used in Fig. |2la) ^zj. (b) The 
largest m:n tongues observed for the chemical conditions used in 
Fig.Hfb) 1 29]. The homogeneous frequency in both cases is /o = 
0.020 Hz while the spiral frequency is (a) 0.072 Hz, (b) 0.020 Hz. 
Each symbol type represents a different m:n response. The patterns 
(points) within the solid curves respond sub-harmonically with the 
forcing frequency. The bottom plot is taken from l33ll . 



transformed. The complex amplitudes of the filtered system 
give information about the phase distribution of the pattern. 
Qualitatively different patterns were found to exhibit different 
shapes in the complex phase-plane representations |34] [see 
Fig.Ee-h)]. 



IV. PATTERN FORMATION 

We now describe the asymptotic patterns observed (after 
the decay of transients) in the m:l tongues for different forc- 
ing frequencies and amplitudes. The patterns can be divided 
into two categories: those which are resonant with the forc- 
ing and those which are near-resonant. The resonant patterns 



are standing waves which lock to the forcing frequency and 
show m peaks in the phase response. The near-resonant pat- 
terns are traveling waves and spiral waves which do not lock 
to the forcing but still show m peaks in the phase response. 

We differentiate the two types of response using their power 
spectra. If the system is resonant, the response frequency will 
adjust to be a rational ratio of the forcing frequency. For near- 
resonant response, however, the frequency does not adjust to 
be a rational ratio of the forcing frequency. Figure |5l shows 
the power spectra and corresponding histograms of phase an- 
gle for 2: 1 near-resonant patterns (spirals) at three forcing fre- 
quencies near exact resonance with the spiral wave frequency. 
The peak of the subharmonic response does not adjust to // /2 
(vertical line) as the frequency is varied. The histograms of 
the phase, however, indicate that there is a two-phase response 
to the forcing even though there is not exact resonance. In 
contrast. Fig. |5] shows the power spectra of resonant patterns 
(standing waves) for three forcing frequencies near exact res- 
onance with the uniform oscillation frequency. In this case the 
patterns lock to // /2 (shown by the vertical line) even when 
the forcing is detuned from exact resonance. 

We now discuss the different types of m: 1 patterns that we 
observed. 



A. 1:1 and 2:1 patterns 

In the 1:1 region we observe a resonant response. In this 
case the entire pattern of chemical concentration oscillates 
uniformly in space with the forcing frequency, as measured 
for a range of // and / values centered at /q. The shape of the 
1:1 tongue in Fig.|3^a) is different than the shape of the other 
tongues. The other tongues bend toward higher frequency at 
low forcing intensities. Instead we find 1 : 1 uniform patterns at 
frequencies near /o even at very low forcing intensities We do 
find spiral patterns at slightly higher frequencies near the bot- 
tom of the tongue, but we cannot distinguish 1 : 1 spiral waves 
from unforced spiral waves. 

Unlike the 1 : 1 resonant response, for which we observed 
only a single qualitative pattern, several qualitatively different 
patterns were observed inside of the 2: 1 region. In this region, 
the oscillation phase responds to either the first or the second 
forcing cycle, which occurs within a single oscillation cycle of 
the pattern. The 2: 1 patterns are therefore formed from spatial 
arrangements of regions oscillating at the same frequency but 
which differ in phase by K. A description of the different 2:1 
patterns observed in the BZ system and in a forced reaction- 
diffusion model with Brusselator kinetics was given in i2lll . 

B. 3:1 patterns 

In the 3:1 region we observe two qualitatively different 
types of patterns. At low forcing the 3:1 patterns are rotat- 
ing spirals, such as those shown in Fig. At low forcing 
intensity, the spirals have a fairly evenly distributed phase an- 
gle and a nearly circular shape in the complex phase plane. At 
higher forcing intensity, the phase becomes more concentrated 



5 





■ 1 



-7r/2 7r/2 Jt 
phase angle 



FIG. 5: Power spectra and phase histograms for resonant 2:1 
standing-wave patterns. Frames (a)-(c) show the power spectra of 
the response as the forcing frequency is varied across the tongue. The 
power spectra, normalized to the forcing frequency ff, show that the 
largest response is at exactly half the forcing frequency (//// = 0.5), 
as indicated by the vertical line. Frames (d)-(f) show histograms of 
the phase of the pattern near the peak response. The histograms 
have two peaks corresponding to concentrations of the pattern in 
regions that are separated in phase by n. Forcing frequency: (a,d) 
// = 0.0333 Hz, (b,e) // = 0.0357 Hz, (c,f) // = 0.0416 Hz. 



in three phases, and the shape in the complex phase plane be- 
comes more triangular. This trend is observable in Fig.0b,c). 

An abrupt transition from traveling wave patterns to what 
appear to be standing wave patterns is observed in the 3:1 res- 
onance region as the forcing amplitude was increased. The 
transition between these pattern types was observed at a fixed 
forcing frequency of 0.075 Hz as / was increased past roughly 
460 W/m2, and was also observed for fixed forcing amplitudes 
in a range of 300 to 400 W/m2 when // was increased past 
roughly 0.065 Hz. The experimental resolution is not enough 
to determine the functional form of the transition. 

The 3:1 standing wave patterns with stationary or nearly 
stationary fronts consist of irregularly shaped domains dif- 
fering in phase by 27i/3 (see Fig. 0. Often the fronts are 
rough, i.e., have short wavelength modulations that appear sta- 
ble over a hundred oscillation cycles of the pattern, as can be 
seen in the standing wave pattern pictured in Fig.[nc). 

The fronts in the standing-wave patterns are either station- 
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FIG. 6: Power spectra and phase histograms for near-resonant 2:1 
spiral wave patterns. Frames (a)-(c) show power spectra of the re- 
sponse as the forcing frequency is varied across the tongue. The 
power spectra, normalized to the forcing frequency ff, show that 
the largest response is nearly half (but not exactly half) the forc- 
ing frequency when the forcing is detuned from exact resonance 
(//// = 0.5 as indicated by the vertical line); thus the spiral wave 
pattern does not lock to the forcing frequency. The peaks at //// = 
1 in (a) and (c) include power from the forcing light reflected from 
the reactor face into the camera. Frames (d)-(f) show histograms 
of the phase of the pattern near the peak response. The histograms 
have two peaks which indicates that the pattern is mostly concen- 
trated near one of two phases; the two phases are separated by it. 
Forcing frequency: (a,d) // = 0.0500 Hz, (b,e) // = 0.0556 Hz, (c,f) 
// = 0.0625 Hz. 



ary or propagate on a time scale orders of magnitude larger 
than the unforced spiral period. If the latter is the case, the pat- 
terns are not precisely standing waves, and over a longer time 
could evolve to other patterns such as large and slowly rotating 
three-phase spiral waves. If so, these patterns could be under- 
stood to be observations of the 3:1 patterns predicted by the 
forced complex Ginzburg-Landau (CGL) equation 113. Bal . 



C. Other patterns 

Rotating four-phase spiral patterns, e.g. Fig. ^d,h), are 
the only pattern type observed in the 4: 1 resonance region of 
the BZ experiments. A detailed description of 4:1 resonance 
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FIG. 7: Patterns observed in the 3:1 resonance region for dif- 
ferent forcing amplitudes (the patterns are shown after filtering). 
The response is shown in the x — y plane (left) and complex phase 
plane (right). (a,b) Spiral waves in the bottom of the 3:1 region, 
/ = 269W/m^, // = 0.1251 Hz. (c,d) Spiral waves in the middle 
of the 3:1 region, / = 382W/m^, // = 0.0667 Hz. (d,e) Standing 
wave patterns found at the top of the 3:1 region, / = 863W/m^, 
// = 0.0769 Hz. 



in the forced CGL equation and in the FitzHugh-Nagumo 
and Brusselator reaction-diffusion models was given in i34ll . 
In those models a bifurcation from four-phase traveling pat- 
terns to two-phase 4:1 resonant standing-wave patterns was 
found 1 13]. This bifurcation has not been observed in our ex- 
periments, perhaps because the forcing light intensity / avail- 
able may be insufficient to reach the bifurcation. Another pos- 
sibility is that the patterns observed were at a frequency near 
four times the spiral wave frequency /, instead of near four 
times the uniform oscillation frequency /q. In the experimen- 
tally accessible range of /, we have observed only 4: 1 spiral 
wave patterns and none of the more complicated pattern be- 
havior found in the 4:1 forced CGL model ll36ll . 

Patterns such as the 2: 1 Bloch fronts and spiral waves dis- 
cussed in 1 35^ and the 4:1 rotating spirals discussed 
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FIG. 8: BZ spirals observed in the 5:1 resonance region at / = 661 
W/m^, ff = 0.100 Hz. The patterns are filtered to keep only the 
frequencies shown by the gray band in (c). (a) Phase patterns in 
X — v plane, (b) Phase in complex plane, (c) Average temporal power 
spectrum, (d) Histogram of phase angle. 

in f]3\ are not resonant since they are all traveling patterns. 

Finally, we present examples of near resonant 5:1 and 6:1 
patterns in Fig.|8land Fig.|9l respectively. 

V. REACTION-DIFFUSION MODEL 

As a model for a periodically forced oscillatory system 
we use a version of the FitzHugh-Nagumo (FHN) reaction- 
diffusion equations 

Uf — M — M^— V + V^M, (la) 

V; = e[M- (fli +Ysinco/f)v] +5V^v, (lb) 

where the fields u(x^y) and v(x^y) represent concentrations 
of chemicals in a simple model of a chemical system. We 
add explicit time dependence to the FHN system as parametric 
sinusoidal forcing with amplitude y and frequency co/. The 
parameter e is the ratio of the time scales of u and v and 5 
is the ratio of the diffusion rates of u and v. In the following 
we fix the parameters £ = 0.1,5 = 0.1,fli ~ 0.5, and vary the 
forcing frequency and amplitude. 

In the absence of forcing, y = 0, the equations have a spa- 
tially uniform solution m = y = 0. The parameter e con- 
trols the stability of this solution. When e>2, u — v — Q 
is stable, and at e = 2, there is a Hopf bifurcation to uni- 
form oscillations. Beyond the Hopf bifurcation Eqs. Q also 
support traveling phase waves. Our numerical investigations 
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FIG. 9: BZ spirals observed in the 6: 1 resonance region at / = 475 
W/m^, ff = 0.1001 Hz. The patterns are filtered to keep only the 
frequencies shown by the gray band in (c). (a) Phase patterns in the 
x — y plane, (b) Phase in complex plane, (c) Average temporal power 
spectrum, (d) Histogram of phase angle. 



FIG. 10: Spiral wave in Eq. Q with periodic forcing at nearly twice 
the spiral wave frequency. The patterns are filtered to keep only the 
frequencies shown by the gray band in (c). (a) Phase pattern in x — 
y plane, (b) Phase in complex plane, (c) Average temporal power 
spectrum, (d) Histogram of phase angle. Parameters: y= 0.5, CO/ = 
0.474, x= [0,256], >-= [0,256]. 



are conducted in the parameter range where uniform oscilla- 
tions and phase waves both exist. In two space dimensions 
phase waves typically form into rotating spirals, each one or- 
ganized around a core where the amplitude of oscillations is 
zero. For the parameters above, the spiral wave frequency 
(cOi « 0.237) is faster than the homogeneous oscillation fre- 
quency (cOi « 0.215); once formed, spiral waves spread to fill 
the entire system. 



A. Periodic forcing and data analysis 

A sinusoidal parametric forcing, homogeneous in space, is 
applied by choosing a nonzero y parameter in Eq. Q. As 
in the BZ experiment, when the forcing amplitude is high 
enough the system can lock at rational ratios of the forcing 
frequency. Frequency locking of spatially uniform solutions 
to Eq. Q occurs in tongue-shaped regions in the (Of —y pa- 
rameter plane. A complete diagram of the pattern-forming 
tongues was not computed for the FHN equation with this 
forcing scheme. We have only investigated certain resonances 
to compare with the BZ chemical experiment. Our numerical 
investigations show that the size and shape of different m:n 
tongues depends on the exact form of the parametric forcing 
in Eqs. Q f3^. Tongue diagrams obtained by forcing other 
terms of the FHN model were presented in studies of locking 
to uniform oscillations in the oscillatory FHN model |4d,'4l||. 
A diagram for the 2: 1 tongue of a periodically forced Brusse- 



lator reaction-diffusion system was given in Ref. 12 IB . 

The data from the numerical solutions of the forced FHN 
equation are processed to extract phase information, as in the 
experimental system. Every 1 .4 time units, we store the values 
of u{x,y) and v{x,y) at each computational grid point Xi,yj. 
The discrete Fourier transform is applied to the time variable 
of u{x,y,t) to get the frequency response u{x,y,(i)) for each 
point in the pattern. The averaged power spectrum of the sig- 
nal, 



Pi(0) = 



1 



N,Ny 



(2) 



where Nx and A', are the number of grid points in the x and y 
directions, is then examined to determine the system response. 
The response frequency is isolated from the signal u using a 
box filter centered at the response frequency, co, with width A. 
The filtered signal is then inverse Fourier transformed, which 
gives the response in time, a{x,y,t). 

Figure [Tol shows an example of a 2:1 spiral wave with 
the peak frequency response at (X)//2 displayed clearly in the 
power spectra. The filtered signal a inx—y domain is shown 
in Fig. llOr a). and Fig. llOf b) shows the same data plotted in 
the complex phase plane. The width of the filter is shown in 
the power spectrum by the gray band (see Fig. llOf c)). The 
phase plane shows that different parts of the spatial domain 
are in different relative phases, all oscillating at the same fre- 
quency. The phase is not uniformly distributed but has peaks 
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near two phases that become apparent in the histogram of the 
phase angle shown in Fig. llOf d). 



B. Pattern Formation 

The FHN equations ([l) have two intrinsic frequencies, the 
uniform oscillation frequency, coo, and the spiral wave fre- 
quency cOj. For some choices of parameters (e.g. 5=1) these 
two frequencies are the same, but for the parameters chosen 
in this study the two frequencies differ (cOj > COo). Because of 
this there are two possible different resonant response condi- 
tions: when the forcing is a rational multiple of either coq or 
G)s- In the following we will show how the FHN equations Q 
respond to forcing in both of those cases. In the BZ experi- 
ment this distinction is harder to make. 



C. Patterns at m:n response of the uniform oscillation 
frequency 

When the forcing is a rational multiple of the uniform os- 
cillation frequency cOq, spatially uniform solutions of Eq. Q 
are found for a range of forcing frequency Oif and amplitude 
y. The m:n frequency-locked solutions form tongue shaped 
regions in the coy — y parameter plane (Arnol'd tongues). The 
shape of the resonant tongues depends on the exact form of 
the forcing |39]. 

As in the BZ experiment, patterns may form in the 
frequency-locked tongues. Resonant pattern solutions con- 
sist of standing waves connecting regions of different phases. 
For example, in the 2:1 resonance, standing waves consist of 
fronts between regions in space that are oscillating at the same 
frequency but out of phase by K. The fronts must be station- 
ary for the pattern to be strictly frequency locked, since any 
motion indicates that the phase is drifting and thus, at least in 
the vicinity of a front, the frequency is also slowly changing. 

A 2:1 resonance is found when the forcing frequency is 
nearly twice the uniform oscillation frequency, cOj ~ 2a)o. 
For sufficiently high forcing amplitude, the system frequency 
locks at (Of/2, in one of two phases separated by K. Pat- 
terns form from the two phases. At low forcing amplitude 
the pattern is a two phase rotating spiral wave and thus is 
near-resonant but not strictly frequency-locked since the spi- 
ral rotates slowly (relative to the forcing period). At higher 
forcing amplitude the pattern is standing waves. The spiral 
waves and standing waves are similar to those found in the 
BZ experiment and in forced complex Ginzburg-Landau and 
Brusselator models 1 14, 21, 22, 42|. 

In the 3:1 resonance the system responds at one-third the 
forcing frequency, (Or = K>//3, and the patterns consist of 
spatial regions of three locked phases. In our exploration in 
forced FHN model we find that the three locked phases orga- 
nize into rotating spiral waves. 

Patterns in a forced 4: 1 FHN model consist of four-phase 
spiral waves and two-phase standing waves and are discussed 
in detail in Ref. fll . 



The domain in parameter space where frequency locked 
patterns (standing waves) exist is different from that of fre- 
quency locked uniform solutions. Recently it was discovered 
that frequency locked standing wave patterns can exist out- 
side of the resonant tongue of spatially uniform solutions or 
that s pati al instabihties can reduce the range of resonant pat- 
terns glH- 

D. Patterns at m:n response of the spiral wave frequency 

When the forcing frequency is close to a rational multiple of 
the spiral wave frequency, the spiral does not frequency lock 
to the forcing but still shows a near-resonant response. This 
can be seen in the nonuniform phase distribution of the forced 
spiral waves, with m peaks for waves forced near an m:n reso- 
nance. For example, forcing the spiral wave at approximately 
twice the spiral frequency, CO/ ~ 2(0,, causes the spiral to re- 
spond by shifting the relative oscillation phases within the 
spiral to be concentrated near two phases, as shown in Fig- 
ures[l3d)-(f). 

Figure [21 shows the spiral wave response when the forc- 
ing frequency is scanned through the spiral frequency. The 
maximum response in the power spectrum is not (Of/2 when 
the forcing frequency is not exactly twice the spiral wave fre- 
quency; the pattern is not frequency locked but it is near- 
resonant. In contrast to a quasi-periodic response farther away 
from resonance i2lll . near-resonant patterns exhibit a nonuni- 
form distribution in the histograms of the phase angle. This 
distribution is farthest from uniform when the forcing is clos- 
est to (Of/2, see Fig.Od)-(f). 

Forcing near other resonances of the spiral wave frequency 
also results in a near-resonant response, with the number of 
peaks in the phase histogram corresponding to the m:n reso- 
nance. In the 3:1 resonance we observe three peaks and in 
the 4: 1 resonance we find four peaks. Other spiral resonances 
such as 5: 1 and 6: 1 can also be found. These results are similar 
to observations near-resonant spirals in the of the BZ experi- 
ment. 

To characterize the effect of the forcing on the spiral wave 
pattern, we measured the deviation of the phase from a uni- 
form distribution. For an unforced spiral wave, the histogram 
of the phase near the spiral frequency is flat, indicating that 
the phase is uniformly distributed between —K and K. When 
the spiral wave is nearly 2:1 resonant, the histogram shows 
two peaks that are separated by K in the phase distribution. 
Figures [T2l d)-(f) show. 

The nonuniform phase response was measured by the chi- 
square statistic relative to the uniform distribution. The phase 
data at k = 256 x 256 computational grid points was binned 
into 100 equal size bins between — 7t and K. The chi-square 
statistic is 

X^^E^. (3) 

where A^,- is the value in bin / and E = k/lQQ is the expected 
value. Figure shows the dependence of on the forcing 
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FIG. 1 1 : Resonant standing-wave response in the 2: 1 forced FHN 
Eq. Q. (a)-(c) Power spectra for tliree different forcing frequencies 
near 2a)o. Tiie frequencies are normalized to tfie forcing frequency 
coy. Tlie peak subliarmonic response is exactly at half the forcing 
frequency, even when the forcing is not exactly 2coo! as indicated by 
the vertical line at co/coj = 0.5. (d)-(e) Histograms of the distribution 
of phase angle for the spiral waves in (a)-(c). The two peaks in the 
distribution show that the pattern response is primarily in two phases 
separated by an angle of K. Parameters: y = 3.0, (a) cOy = 0.349, (b) 
CO/ = 0.370, (c) CO/ = 0.419. 

amplitude y when the forcing frequency is near 2cl)s. The 
value increases exponentially as the forcing amplitude is in- 
creased from near zero. The nonzero value of at y = is 
due to small fluctuations in the phase of our finite size sample. 



VI. DISCUSSION 

The BZ chemical system in an open gel-reactor was used to 
identify multiple tongues, each with a different m:n resonance, 
in the forcing frequency-amplitude parameter plane. Such a 
phase diagram has not been previously reported for a spatially 
extended oscillatory system. The resonance tongues are found 
to be ordered in a Farey sequence, similar to the Devil's stair- 
case ordering of resonance tongues for two coupled oscilla- 
tors 1 30 1 and for the homogeneous BZ reaction ll3lll32ll . 

The diffusively coupled oscillations we measure respond to 
external forcing either resonantly or at near resonance (quasi- 



FIG. 12: Response of a 2:1 forced spiral wave in the FHN Eq. G). 
(a)-(c) Power spectra for three different forcing frequencies (a) 
C0//2 < cOj, (b) C0//2 ^ cOj, (c) C0//2 > cOj. The solid line indi- 
cates (Of/2. When the forcing is not exactly in resonance the peak 
response in the power spectra differs from C0//2 indicating that the 
spiral is not frequency locked, (d)-(e) Histograms of the distribu- 
tion of phase angle for the spiral waves in (a)-(c). The two peaks 
in the distribution show that the spiral wave responds to the forcing 
by redistributing the internal relative distribution of the phase even 
though the pattern is not frequency locked. Parameters: Y= 0.3, (a) 
CO/ = 0.465, (b) CO/ = 0.474, (c) CO/ = 0.493. 



periodically but with an m-peaked phase distribution). The 
resonant patterns are standing waves that frequency lock to a 
m:n ratio of the forcing frequency. In this case, a power spec- 
trum of the resulting resonant pattern shows a single primary 
peak at // / m, along with its higher harmonics, and the phase 
distribution has m peaks shifted by 2n/m. The near-resonant 
patterns are traveling waves which do not lock to a ratio of the 
forcing frequency but have a response near // /m. However, 
the phase distribution still shows m peaks. This near-resonant 
quasi-periodic behavior is different from quasi-periodicity far- 
ther away from resonance, where patterns have a flat phase 
distribution i2lll43l . 

Both the resonant and near-resonant behavior are also ob- 
served in a FitzHugh-Nagumo reaction-diffusion model with 
sinusoidal periodic forcing, similar to the experiments. 
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